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Abstract 

In this paper we investigate the elastic properties of the vortex lattice for 

a superconducting film of finite thickness. We derive an analytic expression 

for the compression modulus. The shear modulus is evaluated numerically by 

using both the Pearl interaction potential, valid in the limit of very thin film, 

and a potential for films of arbitrary thickness. A comparative study of the 

shear moduli is carried out. 
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The problem of a vortex emerging perpendicularly to a surface of a superconductor was 
first considered by Pearl.B'i He pointed out that the vortex- vortex interaction potential at the 
surface of superconductors follows the power law 1/r (1/k in Fourier space) at large distances. 
The Pearl interaction potential has been widely used to investigate the equilibrium!! and the 
elastic^ properties of the vortex lattice in superconducting films. 

In this work we will show that in fact the Pearl interaction potential is capable of describ- 
ing satisfactorily the superconducting properties of a very thin superconductor. However, 
we will show that as the thickness of the film approaches the London penetration depth A, 
the use of this potential may underestimate significantly the value of the shear modulus for 
sufficiently low induction. 

It has been shown elsewhere! that the energy of an ensemble of interacting vortices inside 
a superconducting film of finite thickness d and the energy of the stray fields in the vacuum 
is given by 
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where k = (k x ,k y ), k = \Jk 2 + k 2 , a = \/k 2 + A -2 , and $o is the quantum flux. The 
structure factor is given by 

S(k) = ]T e ikR ' . (2) 

i 

Here Rj is the position of the zth vortex line. Note that this result should be valid for 
an ensemble of distorted vortices, that is, the positions of the vortices do not necessarily 
correspond to the equilibrium configuration. The first term inside Eq. (H) represents the 
interaction energy of the vortex lines as if the surfaces were absent. The second term repre- 
sents the surface energy associated to the magnetic energy of the stray field at the vacuum. 
Note that for k small (large r), a 2 ~ 1/A 2 . Thus, the surface energy goes as $Q/87r 2 r. 
This is the Pearl result for vortices emerging from a semi-infinite isotropic superconductor!] 
Another interesting particular case of Eq. ([!]) is the limit of a very thin film, <i — *> 0, and k 
small. In this limit, from Eq. (J!]) it is straightforward to show that 

where E = ($ /47rA) 2 , A = 2X 2 /d. This is precisely. the energy of an ensemble of interacting 
vortices in a very thin film first obtained by Pearl.@ However, his derivation is supposed to 
be valid for any k. In this work we use Eq. (^) with no restriction on k. 

The integrand of Eq. ([I]) contains two terms. The bulk term follows the power law 1/k 2 
for large k. The integrand of Eq. (|3|) has the same behavior as its corresponding one. As 
a result, the self-energy contributions to both Eqs. (HD and @ diverge. The divergence is 
logarithmic. This is so because London theory neglects the size of the vortex-core. To remove 
this short-length scale divergence we use the Gaussian cutoff. This regularization procedure 
consists in multiplying the integrand of Eq. (JiJ) and (^) by a factor e~ 2 ^ k , where £ is the 
coherence length.0 In the first case, only the bulk term needs a cutoff.! The London picture 
is valid in the limit of low induction and very strong type II superconductor, k — A/£ ^> 1. 
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Let us suppose that initially the vortex lines are located at the equilibrium positions 
R° = (Xf,Y®). Let us denote by Uj = u(R°) the displacement of the ith vortex line from 
its equilibrium position. The new position of the ith vortex line is then given by Rj = R°+u,. 
The elastic energy of the distorted vortex lattice is most conveniently written in k-space by 
expressing the displacement in terms of its Fourier transform 

L <SH u(k)ekR; ' • (4) 

where the integration is taken over the first Brillouin zone (BZ) of the vortex lattice. Upon 
taking this into Eq. ([I]) and expanding it up to second order in Uj we obtain for the excess 
free energy 

1 r dk 2 1 

6F = 2 J J2nj ^Mk)$ a/3 (k)M-k) , (5) 

with (cn,/3) = (x,y). Here B = §o/A is the induction, A is the area of a unit cell of the 
vortex lattice. The factor l/d— j_ n d / d dk z / (2tt) in the integrals was introduced deliberately 
to resemble the three dimensional corresponding expressions. The coefficients $ Q/ g(k), called 
elasticity matrix, are real, symmetric, and periodic in k-space. One has 

$ Q/3 (k) = ^ £ [(k + Q) a (k + Q)^(k + Q) - Q a Q p V(Q)\ , (6) 

47r Q 

where Q are the reciprocal lattice vectors. Although the last equation is valid for any type 
of vortex lattice, we will use a triangular lattice. The basis vectors of the reciprocal lattice 
are given by 

where a 2 = 2$ /v / 3-B. Then, the vortex positions in reciprocal space are given by Q = 
Qmn = ^Qi + m Q2, with m, n integers. 

The interaction potential V(k) for a film of arbitrary thickness is 



V(k) 



A 2 a 2 

and in the limit of Pearl is 



1 + 



d\ 2 ka(k + a) 



(8) 



y W = uT^. ■ P) 

Since we have assumed that the vortex lines are straight and parallel to each other, 
there will be only two elastic constants connected to the excitation of the lattice. One is 
the compression modulus and the other one is the shear modulus. Usually, shear is softer 
than compression. Let us first neglect the discreteness of the lattice by considering only the 
Q = contribution to Eq. (^). Within this continuum approximation, only the compression 
modulus is obtained. One has 
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c n (k) = |V(k) . (10) 

Note that in the long wave length limit, k — > 0, the interaction potential diverges, both 
for the full expression and within the Pearl limit [see Eqs. (f|) and @]. As a result, in this 
local limit the compression modulus diverges. Thus, as the vortex density becomes very 
small, the energy cost to compress the lattice is very large. 

The evaluation of the more important shear modulus is possible only if one goes beyond 
the continuum limit. Over the first Brillouin zone the shear modulus is nearly a constant, 
that is, it does not contain any significant non-locality like the compression modulus.i Thus, 
either of the following definitions is sufficient to describe the shear deformations 

c 66 = hm — , c 66 = hm — . (11) 

fc->o k 2 k^o k 2 

In what follows, the induction B is in units of the upper critical field H c2 , b = B/H c2 \ 
lengths are in units of A; the Ginzburg-Landau parameter used is k = 50. By using these 
parameters we evaluated numerically Cg 6 and c^q, where the superscript p means that the 
Pearl potential was used to calculate the shear modulus, whereas with no superscript we 
used the full expression of Eq. (H). Fig. [I] shows the ratio (%q/c§§ as a function of the 
film thickness. As one can see from this figure, for sufficiently low induction the difference 
between both shear modulus grows as d increases. Note that this occurs for values slightly 
larger than £. However, these differences tend to disappear as the induction increases. 

In Fig. ^|we also show a plot of the ratio Cg 6 /c 66 as a function of the film thickness, where 
the superscript b stands for bulk. Note that as the induction increases the shear modulus 
the bulk shear modulus becomes approximately equal to full shear modulus. This was first 
experimentally observed by Fiory.tl 

The conclusion we draw from this scenario is that the Pearl potential is indeed valid in 
the limit of very thin film. However, the study of the elastic properties of the vortex lattice 
in superconducting films, of thickness slightly larger than £ and the induction very close to 
the lower critical field H c \, might become non-reliable by employing this potential. In this 
situation it should be more convenient to use the full expression for the interaction potential, 
but for films of thickness not much larger than A, because in this limit bending of the vortex 
lines may become important. 
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FIGURES 



FIG. 1. The ratio c^q/cqq for several values of the induction, (a) b = 0.0005 (continuous line); 
(b) b = 0.005 (dot-dashed line); (c) b = 0.01 (dashed line). 

FIG. 2. The ratio c b m /cQQ for several values of the induction, (a) b = 0.001 (continuous line); 
(b)b = 0.002 (dot-dashed line); (c) b = 0.005 (dashed line). 
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